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A. PDFs and Nuclear Corrections 

Parton distribution functions (PDFs) are of supreme 
importance in contemporary high energy physics as they 
are needed for the computation of reactions involving 
hadrons based on QCD factorization theorems pHlj]. F° r 
this reason various groups prese nt glob al analyses of 
PDFs for protons [J-[ll| and nuclei 12t13| which are reg- 
ularly updated in order to meet the increasing demand 
for precision. The PDFs are non-perturbative objects 
which must be determined by experimental input. To 
fully constrain the cc-dependence and flavor-dependence 
of the PDFs requires large data sets from different pro- 
cesses which typically include Deeply Inelastic Scattering 
(DIS), Drell-Yan (DY), and jet production. 

While some of this data is extracted from free protons, 
much is taken from a variety of nuclear targets. Because 
the neutrino cross section is so small, to obtain suffi- 
cient statistics for the neutrino-nuclear DIS processes it 
is necessary to use massive targets (e.g., iron, lead, etc.). 
Therefore, nuclear corrections are required if we are to 
include the heavy target data into the global analysis of 
proton PDFs. 

The heavy target neutrino DIS data plays an impor- 
tant role in extracting the separate flavor components of 
the PDFs. In particular, this data set gives the most 
precise information on the strange quark PDF. As the 
strange quark uncertainty may limit the precision of par- 
ticular Large Hadron Collider (LHC) W and Z measure- 
ments, the nuclear corrections and their uncertainties will 
have a broad impact on a comprehensive understanding 
of current and future data sets. 
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Figure 1: Nuclear correction ratio, F-f^/F^P, as a function 
of x. The parameterized curve is compared to SLAC and 
BCDMS data pHU. 



B. Nuclear Corrections in the Literature 

In previous PDF analyses (25l . [2f| , a fixed nuclear cor- 
rection was applied to "convert" the data from a heavy 
target to a proton. As such, these nuclear correction fac- 
tors were frozen at a fixed value. They did not adjust for 
the Q 2 scale or the physical observable (F2, F3, 3^), 
and they did not enter the PDF uncertainty analysis. 

While this approach may have been acceptable in the 
past given the large uncertainties, improvements in both 
data and theory precision demand comparable improve- 
ments in the treatment of the nuclear corrections. 

Figure [T] displays the F^ e /F® structure function ratio 
as measured by the SLAC and BCDMS collaborations. 
The SLAC/NMC curve is the result of an A- independent 
parametrization fit to calcium and iron charged-lepton 
DIS data [IB-H, H3 ■ This parameterization was used to 
"convert" heavy target data to proton data, which then 
would be input into the global proton PDF fit. 1 The 
SLAC/NMC parmeterization was then applied to both 
charged-lepton-nucleus and neutrino-nucleus data, and 
this correction was taken to be independent of the scale 
Q and the specific observable {F2, F3, ...}. Recent work 
demonstrates that the parameterized approximation of 
Fig. [1] is not sufficient and it is necessary to account for 
these details (HHIJ. 



C. Outline 

In this paper, we present a new framework for a global 
analysis of nuclear PDFs (NPDFs) at Next-to-Leading- 
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Figure 2: We display the ^4-dependent coefficients Ck(A), 
k — {1, 5}, for the up-valence (top) and down-valence PDF 
(bottom) as a function of the nuclear A. The dependence of 
the coefficients Ck(A) is shown by the following lines: ci- solid 
(red) line, C2- long dashed (blue) line, C3- dashed (green) line, 
C4- dash-dotted (magenta) line, c 5 - dotted (brown) line. 

Order (NLO). An important and appealing feature of 
this framework is that it naturally extends the proton 
analysis by endowing the free fit parameters with a de- 
pendence on the atomic number A. This will allow us to 
study proton and nuclear PDFs simultaneously such that 
nuclear correction factors needed for the proton analysis 
can be computed dynamically. 

In Section [TTl we outline our method for the analysis, 
specify the DIS and DY data sets, and present the \ 2 of 
our fit. In Section Hill we compute the nuclear correction 
factors {F( e /F®) for the fit to the l±A and DY data. 
In Section IIV1 we compare these results to the nuclear 
correction factors (Ff e /Ff) from the vA fit of Ref. (30j . 
Finally, we summarize our results in Section Ivl 

II. NPDF GLOBAL ANALYSIS FRAMEWORK 



Technically, the heavy target data was scaled to a deuteron tar- 
get, and then isospin symmetry relations were used to obtain 
the corresponding proton data. Deuteron corrections were used 
in certain cases. 



A. PDF analysis framework 

In this section, we present the global analysis of NPDFs 
using charged-lepton DIS (I A) and Drell-Yan data to 
extend the analysis of Ref. [27j for a variety of nuclear 
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Figure 3: We display the a) xu(x) and b) xd(x) PDFs for a selection of nuclear A values ranging from A — {1,207}. We 
choose Qo — 1.3 GeV. The different curves depict the PDFs of nuclei with the following atomic numbers (from top to bottom 
at x = 0.01) A = 1,2, 4, 8, 20, 54, and 207. 



targets. This analysis is performed in close analogy with model and the \ 2 analyses as outlined in Ref. |30|. The 
what is done for the A = 1 free proton case |45||. We input distributions are parameterized as 
will use the general features of the QCD-improved parton 



xf k (x,Qo) = c Q x Cl (l-x) C2 e C3X {l + e Ci x) C5 k = u v ,d v ,g,u + d,s,s , (1) 

d(x,Q )/u{x,Q ) = c .t Ci (1-z) C2 + (1 + c 3 x)(l - x) c * , 



at the scale Qo = 1.3 GeV. Here, the u v and d v are 
the up- and down-quark valence distributions, u, d, s, s 
are the anti-up, anti-down, strange and anti-strange sea 
distributions, and g is the gluon. 

In order to accommodate different nuclear target ma- 
terials, we introduce a nuclear v4-dependence in the Ck 
coefficients: 

c k -> c k (A) = c M + c M (1 - A- Ck - 2 ) , k = {1, . . . , 5} . 

(2) 

This ansatz has the advantage that in the limit A — > 1 we 
have Ck(A) —> Ckfl', hence, Ck,o is simply the correspond- 
ing coefficient of the free proton. Thus, we can relate the 
Cfc,o parameters to the analogous quantities from proton 
PDF studies. 

It is noteworthy that the x-dependence of our input 
distributions f^ A (x,Qo) is the same for all nuclei A; 
hence, this approach treats the NPDFs and the proton 
PDFs on the same footing. 2 Additionally, this method 
facilitates the interpretation of the fit at the parameter 
level by allowing us to study the Ck(A) coefficients as 



2 The nuclear analogue of the scaling variable x is defined as x : = 
Ax a where xa = Q 2 /2Pa '<? is the usual Bjorken variable formed 
out of the four-momenta of the nucleus (Pa) and the exchanged 
boson (q), with Q 2 = -q 2 [30ll . 



functions of the nuclear A parameter. 

With this A-generalized set of initial PDFs, we can ap- 
ply the DGLAP evolution equations to obtain the PDFs 
for a bound proton inside a nucleus A, f?^ A [x,Q) . We 
can then construct the PDFs for a general (^4, Z)-nucleus: 

f^M) = \ f? /A (x,Q)+ { -^f: /A (x,Q) (3) 

where we relate the distributions of a bound neutron, 
f^ A (x,Q), to those of a proton by isospin symmetry. 
Similarly, the nuclear structure functions are given by: 

F} A > Z \x,Q) = | F? A {x,Q) + ^1 f:' a (x,Q) . 

(4) 

These structure functions can be computed at next-to- 
leading order as convolutions of the nuclear PDFs with 
the conventional Wilson coefficients, i.e., generically 

F^ z \x,Q)=Y,C ik ®fl A > Z) . (5) 
fc 

To account for heavy quark mass effects, we calculate 
the relevant structure functions in the Aivazis-Collins- 
Olness-Tung (ACOT) scheme 0,113 at NLO QCD 
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F A /F° : 

Observable 


Experiment 


Rcf. 


# data 


D 


NMC-97 


m 


275 


He/D 


SLAC-E139 


m 


18 




NMC-95,re 


m 


16 




Hermes 


m 


92 


Li/D 


NMC-95 


m 


15 


Be/D 


SLAC-E139 


El 


17 


C/D 


EMC-88 


m 


9 




EMC-90 


cm 


2 




SLAC-E139 


El 


7 




NMC-95, re 


m 


16 




NMC-95 


m 


15 




FNAL-E665-95 


m 


4 


N/D 


BCDMS-85 


m 


9 




Hermes 


m 


92 


Al/D 


SLAC-E049 




18 




SLAC-E139 


m 


17 


Ca/D 


EMC-90 


m 


2 




SLAC-E139 


El 


7 




NMC-95,re 


m 


15 




FNAL-E665-95 


m 


4 


Fo/D 


BCDMS-85 


El 


6 




BCDMS-87 


El 


10 




SLAC-E049 


El 


14 




SLAC-E139 


El 


23 




SLAC-E140 


El 


6 


Cu/D 


EMC-88 


m 


9 




EMC-93(addendum) 


M 


10 




EMC-93(chariot) 




9 


Kr/D 


Hermes 


m 


84 


Ag/D 


SLAC-E139 


m 


7 


Sn/D 


EMC-88 


m 


8 


Xc/D 


FNAL-E665-92(em cut) 


m 


4 


Au/D 


SLAC-E139 


m 


18 


Pb/D 


FNAL-E665-95 


m 


4 


Total: 






862 



Observable 


Experiment 


Ref. 


# data 


Be/C 


NMC-96 


El 


15 


Al/C 


NMC-96 


El 


15 


Ca/C 


NMC-95 


El 


20 




NMC-96 


El 


15 


Fe/C 


NMC-95 


El 


15 


Sn/C 


NMC-96 


El 


144 


Pb/C 


NMC-96 


El 


15 


C/Li 


NMC-95 


El 


20 


Ca/Li 


NMC-95 


El 


20 


Total: 






279 



Table II: The DIS F^/F^' data sets used in the fit. The table details 
the specific nuclear targets, references, and the number of data points 
without kincmatical cuts. 



pA / pA' . 
DY / DY • 

Observable 


Experiment 


Ref. 


# data 


C/D 


FNAL-E772-90 


El 


9 


Ca/D 


FNAL-E772-90 


El 


9 


Fe/D 


FNAL-E772-90 


El 


9 


W/D 


FNAL-E772-90 


El 


9 


Fe/Be 


FNAL-E866-99 


El 


28 


W/Be 


FNAL-E866-99 


El 


28 


Total: 






92 



Table I: The DIS F^/F-f data sets used in the fit. The table details 
the specific nuclear targets, references, and the number of data points 
without kincmatical cuts. 



Table III: The Drell-Yan data sets used in the fit. The table details 
the specific nuclear targets, references, and the number of data points 
without kincmatical cuts. 



is listed in Tables HI [TTT and lllll there are 1233 data points 
before kinematical cuts are applied. 

The structure of the fit is analogous to that of Ref. [3Cj ■ 
For the quark masses we take m c = 1.3 GeV and mj, = 
4.5 GeV. To limit effects of higher-twist we choose stan- 
dard kinematic cuts of Q cu t = 2.0 GeV, and W cu t = 
3.5 GeV as they are employed in the CTEQ proton anal- 
yses. 3 There are 708 data points which satisfy these cuts. 
The fit was performed with 32 free parameters which 
gives 676 degrees of freedom (dof). 



B. Inputs to the Global NPDF Fit 

Using the above framework, we can then construct a 
global fit to charged-lepton-nucleus (I A) DIS data and 
Drell-Yan data. To guide our constraints on the Ck,o co- 
efficients, we use the global fit of the proton PDFs based 
upon Ref. [27} • This fit has the advantage that the ex- 
tracted proton PDFs have minimal influence from nuclear 
targets. To provide the A-dependent nuclear informa- 
tion, we use a variety of / ± A DIS data and Drell-Yan 
data. The complete list of nuclear targets and processes 



C. Result of the NPDF Fit 

Performing the global fit to the data, we obtain an 
overall x 2 /dof of 0.946. Individually, we find a x 2 /pt of 
0.919 for the F^/Ff measurements of Table [fl of 0.685 
for the FA/FA' measurements of Table HU and of 1.077 



3 For example, see the CTEQ (Coordinated Theoretical- 
Experimental project on QCD) analysis of Ref. |4"J | which 
presents the CTEQ6 PDF sets. 



5 



for the Drell-Yan measurements of Table Mil The fact 
that we obtain a good fit implies that we have devised 
an efficient parameterization of the underlying physics. 

The output of the fit is the set of Ck,i parameters and 
a set of A-dependent momentum fractions for the gluon 
and the strange quark. Using the Cfc , coefficients we can 
construct the A-dependent Ck(A) functions which deter- 
mine the nuclear PDFs at the initial Qq scale: f^{x, Qo). 
As an example, we display the Ck(A) functions in Fig. [2] 
for the case of the up- valence and down-valence distribu- 
tions. 

Finally, we can use the DGLAP evolution equations to 
evolve to an arbitrary Q to obtain the desired f{ L {x,Q) 
functions. In Fig. [3] we display the up- and down-quark 
PDFs at a scale of Qo = 1-3 GeV as a function of x for a 
variety of nuclear- A values. 



III. e ± A NUCLEAR CORRECTIONS 

Nuclear corrections are the key elements which allow 
us to combine data across different nuclear targets and 
provide maximum information on the proton PDFs. As 
the nuclear target data plays a critical role in differenti- 
ating the separate partonic flavors (especially the strange 
quark) , this data provides the foundation that we will use 
to make predictions at the LHC. 



A. Charged-Lepton (^A) Data 

The present nuclear PDF global analysis provides us 
with a complete set of NPDFs f A (x,Q) with full func- 
tional dependence on {x, Q,A}. Consequently, the tra- 
ditional nuclear correction F^ e /F® does not have to be 
applied as a "frozen" external factor, but can now be- 
come a dynamic part of the fit which can be adjusted to 
accommodate the various data sets. 

Having performed the fit outlined in Sec. [ill we can 
then use the f A (x,Q) to construct the corresponding 
quantity F^ e /F® to find the form that is preferred by the 
data. In order to construct the ratio, we use the expres- 
sion given by Eq. 0] for iron and deuterium. This result 
is displayed in Figure HJ-a) for a scale of Q 2 = 5 GeV 2 , 
and in Figure [5]-a) for a scale of Q 2 = 20 GeV 2 . Com- 
paring these figures, we immediately note that our ratio 
F^ e jF® has non-trivial Q-dependence — as it should. 

Figures 0]-a) and [S}a) also compare our extracted 
F£ e /F? ratio with the (Q-independent) SLAC/NMC 
parameterization of Fig ure \T\ and with the fits from 
Kulagin-Petti (KP) [28|, Hi] and Hirai-Kumano-Nagai 
(HKN07) [IH ■ We observe that in the intermediate range 
(x 6~ [0.07, 0.7]) where the bulk of the SLAC/NMC data 
constrains the parameterization, our computed F^ e jF® 
ratio compares favorably. When comparing the different 
curves, one has to bear in mind the following two points. 
First, all curves in principle have an uncertainty band 



which is not shown. Second, the data points used to ex- 
tract the SLAC/NMC curve are measured at different 
Q 2 whereas our curve is always at a fixed Q 2 = 5 GeV 2 
or Q 2 = 20 GeV 2 . In light of these facts, we con- 
clude that our fit agrees very well with other models 
and parametrizations as well as with the measured data 
points. 

It should be noted that the kinematic cuts we employed 
to avoid higher twist effects effectively exclude all data 
points in the high-.T region above x > 0.7. This is re- 
flected by the fact that our curves in Figs.0J-a) and[S]-a) 
stop at x = 0.7. The high- a; region is beyond the scope 
of this paper and will be subject of a future analysis. 

Thus, we find that data sets used in this fit (FA/F®, 
FA/FA\ and crg^/ag^ ) are compatible with the SLAC, 
BCDMS, and NMC data. Additionally, we can go further 
and use our complete set of NPDFs ff{x 1 Q) to compute 
the appropriate nuclear correction not only for F^' e /F^ ) , 
but for any nuclear target (A) for any Q-value, and for 
any observable. We make use of this property in the 
following section where we compute the corresponding 
quantity for a different nuclear process. 



IV. e ± A AND vA NUCLEAR CORRECTIONS 
A. Nuclear Corrections in vA DIS 

In a previous analysis [30| . we examined the charged 
current (CC) neutrino-nucleus DIS process vA [iX. 
and extracted the F^ e JF® ratio. 4 

These results are displayed in Figures HJ-b) and [SJb) . 
The solid line is the result of the global fit (fit A2), 
and this is compared with the previous SLAC/NMC pa- 
rameterization, as well as fits KP and HKN07. The 
data points displayed come from the NuTeV experiment 
[49L l5Q | . The (yellow) band is an approximation of the 
uncertainty of the fits. 

As observed above, the SLAC/NMC parameterization 
is generally consistent with the results of KP and HKN 
as well as our B fit to £^A and DY data. However, the 
A2 fit of Figures 01-b) and [SJb) does not agree with any 
of these three results. We now examine this in detail. 



B. I A and vA Comparison 

The contrast between the charged- lepton (l^A) case 
and the neutrino (vA) case is striking; while the charged- 
lepton results generally align with the SLAC/NMC, KP, 
and HKN determinations, the neutrino results clearly 



4 While Ref. [jjtj extracted the nuclear PDFs using only the NuTeV 
neutrino— iron DIS data, Ref. |27|| demonstrated that the Chorus 
neutrino— lead DIS data|5ll| was consistent with the NuTeV data 
set. 



(i 




X 



(a) (6) 

Figure 4: The computed nuclear correction ratio, F^ e / ' , as a function of i for Q 2 = 5GeV 2 . Fi gure-a) shows the fit (fit 
B) using char ged - lepton-nucleus (i A) and DY data whereas Figure-b) shows the fit using neutrino- nucleus (vA) data (fit 
A2 from Ref. [301]'). Both fits are compared with the SLAC/NMC parameterization, as well as fits from Kulagin-Petti (KP) 
(Ref. 28, 29]) and Hirai et al. (HKN07), (Ref. fl2l|). The data points displayed in Figure-a) are the same as in Fig.[T]and those 
displayed in Figure-b) come from the NuTeV experiment [49l. |5C||. 
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Figure 5: Same as Figure rjfor Q 2 = 20GeV 2 



yield different behavior in the intermediate x-region. We 
emphasize that both the charged-lepton and neutrino re- 
sults are not a model — they come directly from global 
fits to the data. To emphasize this point, we have su- 
perimposed illustrative data point in Figu res ®b) and 
|SJ-b); these are simply the vA DIS data (43. l50| scaled by 
the appropriate structure function, calculated with the 
proton PDF of Ref. (H- 

The mis-match between the results in charged-lepton 
and neutrino DIS is particularly interesting given that 
there has been a long-standing "tension" between the 
light-target charged-lepton data and the heavy-target 
neutrino data in the historical fits (52l . [53| . This 
study demonstrates that the tension is not only between 
charged-lepton light-target data and neutrino heavy- 
target data, but we now observe this phenomenon in 
comparisons between neutrino and charged-lepton heavy- 
target data. 

There are two possible interpretations of this result. 
1. There is, in fact, a single "compromise" solution for 



the F^ e /F® nuclear correction factor which yields 
a good fit for both the vA and £ ± A data. 

2. The nuclear corrections for the t^A and vA pro- 
cesses are different. 

Considering possibility 1), the "apparent" discrepancy 
observed in Figures 0] and [S] could simply reflect uncer- 
tainties in the extracted nuclear PDFs. The global fit 
framework introduced in this work paves the way for 
a unified analysis of the £^A, DY, and vA data which 
will ultimately answer this question. Having established 
the nuclear correction factors for neutrino and charged- 
lepton processes separately, we can combine these data 
sets (accounting for appropriate systematic and statisti- 
cal errors) to obtain a "compromise" solution. 5 



5 While it is straightforward to obtain a "fit" to the combined neu- 
trino and charged-lepton DIS data sets, determining the appro- 
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If it can be established that a "compromise" solu- 
tion does not exist, then the remaining option is that 
the nuclear corrections in neutrino and charged-lepton 
DIS are different. This idea has previously been dis- 
cussed in the literature [28|, H^, [54|. We note that 
the charged-lepton processes occur (dominantly) via 7- 
exchange, while the neutrino-nucleon processes occur 
via lU^-exchange. Thus, the different nuclear correc- 
tions could simply be a consequence of the differing 
propagation of the intermediate bosons (photon, W) 
through dense nuclear matter. Regardless of whether 
this dilemma is resolved via option 1) or 2), understand- 
ing this puzzle will provide important insights about 
processes involving nuclear targets. Furthermore, a 
deeper understanding could be obtained by a future high- 
statistics, high-energy neutrino experiment using several 
nuclear target materials (5^413 • 

V. CONCLUSIONS 

We presented a new framework to carry out a global 
analysis of NPDFs at next-to- leading order QCD, treat- 
ing proton and nuclear targets on equal footing. Within 
this approach, we have performed a x 2_aiia ly s i s of nu- 
clear PDFs by extending the proton PDF fit of Ref. H3 
to DIS l±A and Drell-Yan data. The result of the fit 
is a set of nuclear PDFs which incorporate not only the 
{x, Q}-dependence, but also the nuclear- A degree of free- 
dom; thus we can accommodate the full range of nuclear 
targets from light (A = I) to heavy (A = 207). We find 
a good fit to the combined data set with a total x 2 /dof 
of 0.946 demonstrating the viability of the framework. 

We have used our results to compute the nuclear cor- 
rections factors, and to compare these with the results 
from the literature. We find good agreement for those 
fits based on a charged-lepton data set. 

Separately, we have compared our nuclear corrections 
(derived with a charged-lepton data set) with those com- 
puted using neutrino DIS (vA — > fiX) data sets. Here, 
we observe substantive differences. 

This fit is novel in several respects. 

• Since we constructed the nuclear PDF fits analo- 
gous to the proton PDF fits, this framework allows 



a meaningful comparison between these two distri- 
butions. 

• The above unified framework integrates the nuclear 
correction factors as a dynamic component of the 
fit. These factors are essential if we want to use 
the heavy target DIS data to constrain the strange 
quark distribution of the proton, for example. 

• This unified analysis of proton and nuclear PDFs 
provides the foundation necessary to simultane- 
ously analyze l^A, DY and, vA data. This will 
ultimately help in determining whether 1) a "com- 
promise" solution exists, or 2) the nuclear correc- 
tions depend on the exchanged boson (e.g., 7/Z or 
W±). 

The compatibility of the charged-lepton l^A and 
neutrino-nucleus vA processes in the global analysis is 
an interesting and important question. The resolution 
of this issue is essential for a complete understanding of 
both the proton and nuclear PDFs. 
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